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Abstract 

We report a first principles, self-consistent, all electron, linear combination 
of Gaussian orbitals (LCGO) calculation of a comprehensive collection of 
magneto-optical properties of nickel and iron based on density functional the- 
ory. Among the many magneto-optical effects, we have studied the equatorial 
Kerr effect for absorption in the optical as well as soft X-ray region, where it 
is called X-ray magnetic linear dichroism (X-MLD). In the optical region the 
effect is of the order of 2% while in the X-ray region it is of the order of 1% 
for the incident angles considered. In addition, the polar Kerr effect. X-ray 
magnetic circular dichroism (X-MCD) and total X-ray absorption at the L2,3 
edges, soft X-ray Faraday effect at the L2,3 edges have also been calculated. 
Our results are in good agreement with experiments and other first principles 

methods that have been used to calculate some of these properties. 
PACS numbers: 71.25.Pi, 78.20.Ls 
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I. INTRODUCTION 



Although magneto-optical properties of magnetic metals have been known for over a hun- 
dred years,§ it is only in the past couple of decades that vigorous interest has been refocussed 
on this subject, partly due to their potential for application in the technology of high density 
data storage. Because of advances in laser and tunable synchrotron sources, a variety of 
different magneto-optical effects including the magneto-optical Kerr effect (MOKE),i the 
Faraday effect ,0 X-ray magnetic circular dichroism for absorption!'! as well as angle resolved 
photoemission,0'0 and X-magnetic linear dichroism in angle resolved photoemission0lli 
have now been extensively studied. The connection between all of these phenomena and the 
electronic structure of the materials in which they are seen, has long been known.lll'lli How- 
ever, only in recent years has there been an effort to perform first principles band structure 
calculations of some of these effects. One electron band calculations for bulk and surfaces, 
alloys and multilayers have been already reported for the polar MOKE,llHll X-MCD and 
X-MLD,&I1 and the X-ray Faraday effect .i 

This paper reports a first principles, self-consistent, all electron calculation of several 
magneto-optical properties of bulk nickel and iron, using the linear combination of Gaussian 
orbitals (LOGO) method. These include the equatorial Kerr effect. X-ray magnetic linear 
dichroism (X-MLD) in absorption at the 2p edges, the soft X-ray Faraday effect at the 2p 
edges, the X-MCD and absorption spectra at the 2p edges as well as the polar MOKE. 
The equatorial Kerr effect defined here is the difference in absorption of p-polarized light 
incident obliquely on the magnetized sample in the equatorial geometry, when the sample 
magnetization is reversed. When observed at core-level edges this effect can be called pho- 
toabsorption X-MLD of a new kind to distinguish it from the first kind observed by holding 
the photon polarization parallel to the x-axis and rotating the magnetization from the x to 
the z axis,ii the second kind obseved due to the rotation of the photon polarization vector 
with the magnetization held fixedH and the one observed in angle resolved photemission. 

Our motivation in performing these calculations has been two-fold. Firstly, since it is 
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known that these effects depend sensitively on the accuracy of the band structure calculation, 
we wanted to demonstrate that our modified tight binding method that includes spin-orbit 
coupling in a very straightforward manner, can produce results at least as accurate as other 
ab-initio electronic structure approaches that apply either the fully relativistic machinery or 
use more complicated computational procedures to avoid using the Kramers Kronig (KK) 
transformation. Our MOKE results on nickel and iron essentially prove this point. Secondly, 
although experiments on the equatorial Kerr effect have been done previously in the optical 
region,ii no experiments on the above-mentioned photoabsorption X-MLD at the 2p edge 
of nickel and iron have been reported. It would thus be interesting to see how well our 
theoretical results for the 2p edge agree with future experiments. A recent experiment on 
the soft X-ray Faraday effect at the 2p edges of iro has also spurred us into calculating 
this effect using the LCGO method. 

We thus begin Section II by outlining in brief the inclusion of spin-orbit coupling to the 
non-relativistic LCGO method, which is slightly different from the previous LCGO work.ii 
We also give a brief description of the fast and efficient KK transformation method we 
have employed in our analysis. In section III we discuss in detail the results of our first 
principles calculation of the elements of the conductivity tensor and their subsequent use in 
determining the magneto-optical properties mentioned above. We shall also point out how 
our method compares with previous theoretical results and with available data. 



II. METHOD OF CALCULATION 

The band structure calculation with spin-orbit coupling included within the LCGO 
method was first done by Wang and Callawayll for nickel and by Singh, Wang and Callaway 
for iron.isi However, the exchange potential in those calculations was of the Xa kind and 
the set of 38 basis functions was also small. Subsequently, the computational procedure was 
revised considerably,^ with the exchange potential being replaced by the more accurate von- 
Barth and Hedin type as parametrized by Rajagopal, Singhal and Kimball,ii and the basis 
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set expanded to a total of 75 functions for 3d transition metals. The method of evaluating 
the Brillouin zone (BZ) integrals was also improved by using the linear analytic tetrahedron 
method,S'S and improvements were made in calculating the K = Fourier component 
of the Coulomb potential.^ This non-relativistic LOGO method was applied to a host of 
different elemental magnetic as well as non-magnetic solids over the past several years and 
a variety of electronic properties such as the Fermi surface, the optical conductivity and the 
Compton profile were calculated&@ which were in fairly good agreement with experiments. 
For the 3d transition metals the basis set thus consisted of 13 s-type, 10 p-type, 5 d-type and 
1 f-type Gaussian orbitals, using published orbitals based on atomic calculations. However, 
relativistic corrections were completely ignored in these improvements. 

With a view to tackling the anisotropic properties of magnetic metals, we have added 
the spin-orbit coupling, Darwin and relativistic mass-velocity terms to the existing code. 
As is well-known, the principal advantage to using a Gaussian basis in a first principles 
tight-binding method is the ability to use analytic expressions for the overlap and several 
Hamiltonian matrix elements. This is true in the case of the spin orbit interaction term 
as well. We have used essentially the same approximations used earlieiii for the spin-orbit 
coupling term i.e., the principal contribution to the spin-orbit matrix elements is only for 
those in the p — p and d — d blocks. We also used the same central-cell approximation 
wherein we retain only spin-orbit matrix elements between those orbitals centered on the 



same atomic site. The difference between our new approach and theirs is that in Ref. ^ the 
spin orbit matrix elements were evaluated in reciprocal space and a Fourier sum had to be 
performed over K vectors to obtain the real space potential. Since the spin-orbit coupling 
is strongest at smaller distances, this necessitated the use of an asymptotic expansion for 
the handling of larger K values to obtain a reasonable value for the K summation. We have 
avoided this completely by evaluating the spin-orbit matrix elements in real space directly. 
In terms of the Gaussian basis, the expresssion for the spin orbit matrix elements of the 
p — p block is. 
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while for those of the d — d block is 



where 



k 

In the above expressions Np and A^^^ are the product of the normalization constants of 
the appropriate Gaussian orbitals for the p — p and d — d block respectively, rih is the number 
of basis functions, a is the sum of the exponents of the Gaussian orbitals, is the number 
of unit cells and Q is the unit cell volume, g{k) is the weight associated with each k point 
in the BZ, lij is the sum of the orbital quantum numbers of the basis functions i and j, 
W = J2k9W is the total weight, and T{z) is the Gamma function. For the charge density 
appearing in the Eq. we used the resulting wavefunctions of the non-relativistic, self- 
consistent calculation. With spin-orbit couphng included, the new generalized eigenvalue 
equation of order 150x150 was diagonalized in l/16th of the BZ at only 219 k points for 
fee nickel and only 125 k points for bcc iron, fewer than the number of points used in the 
previous calculations. 

Using the results of the band structure calculation, the elements of the conductivity 
tensor may be found from the formulas!!. 
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.00 — UJnl{k) + i6 iO + UJnl{k) + iS 

where / goes over the occupied states and n goes over the unoccupied states and II's are the 
k dependent matrix elements of the momentum operator. 

The usual procedure is to evaluate the real part of a^x and the imagninary part of a^y 
by replacing the Lorentzian by a delta function in the limit of 6 going to zero. It is then 



customary to keep S finite to simulate a finite relaxation time. For calculations of magneto- 
optical properties, it is then necessary to perform the KK transformations to obtain the 
imaginary part of a., and the real part of a.^i The KK integrals are known to suffer from 
problems of slow convergence and the necessity of choosing high cutoff values for energy. One 
approachlli evaluates the original Kubo formula directly with lifetime effects as a parameter. 
This method, although accurate, is extremely computationally intensive. We have instead 
performed the KK transformations through the use of two succesive Fast Fourier transforms. 
This method is commonly used in studies of infrared intensities of liquids.^ The basis of this 
method comes from the well known relation that if F{t) is the Fourier transform of f{uj), 

F{t) = ^ / du f{u) (6) 

and h{uS) is the Hilbert transform of /(cu), 

h{u) = -P duj'^^ (7) 

then 

H{t) = -z sgn(t) F{t) (8) 

where H{t) is the Fourier transform of h{uj), 

1 f°° 

H{t) = ^ dujh{uo)e-'^' (9) 
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Preliminary tests on known Hilbert transforms using this method yield results that are 
accurate to better than 1%. The advantages of this method are immense since with sufficient 
number of points, the method is very accurate and fast. Furthermore, it has the flexibility 
of incorporating the lifetime effects very conveniently by simply multiplying the right hand 
side of Eq. by an exponential damping factor, 

H{t) = -i sgn(t) F{t) e^"^!*! (10) 

The disadvantage of using this method is the same as that which arises when trying to 
evaluate the KK integral directly, namely, the need to have function values for frequencies 



more than twice the range of frequencies of interest. This is particularly troublesome for 
functions that don't quickly die down within the energy range of interest. In our problem 
this requires the values of the momentum matrix elements for states up to 30 eV above the 
Fermi level. In the following section we shall show our results of using this method for the 
elements of the conductivity tensor in the optical as well as X-ray region and compare it 
with non-KK transformed results. 

We have applied the results of the foregoing analysis to determine both the polar and 
the equatorial Kerr effects. For the polar Kerr effect the complex Kerr rotation is given by 
the relation,lll 

(11) 



We have calculated the equatorial Kerr effect in the optical as well as the soft X-ray 
region. It is however well known that in the equatorial geometry, the reflection coefficient 
of p-polarized light at oblique incidence depends on the direction of magnetization because 
of its dependence on the off-diagonal component of the conductivity tensor. Thus it is 
evident that a reversal in magnetization should cause a change in the absorbed intensity 
of p-polarized light. Moreover, since such an effect does not occur for s-polarized light the 
equatorial Kerr effect can be observed^! even using unpolarized light. Past calculation Jllii 
for the equatorial Kerr effect in the optical region have used expressions that are correct 
up to first order in K2, the off-diagonal component of the dielectric tensor. To calculate 
the effect in the X-ray region we need the exact expression for the reflection coefficient of 
p-polarized light incident at an angle of incidence 6', 



(12) 



where = 1 — arn^O/n^ and the complex refractive index n \s = ki + k\/ ki where ki 
is the diagonal element of the dielectric tensor. From Eq. (|I2]) the absorption (1 — |rp) can 
be obtained. 
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Lastly, we also calculated the X-ray Faraday rotation for nickel and iron at the 2p edge. 
Here we have used the standard expression for the Faraday rotation 6'^ for a thickness d, 

ujd 



dp = — Re 
2c 



(13) 

where n'"'' = \Jl + 'ima^'^oj) / uj and a^'\(jj) = Cxxi^) ± icrxy{uj). 



III. RESULTS AND DISCUSSION 



A. Elements of the conductivity tensor 



Using the results of the band structure calculation, we calculated clxi^) and o"^y(co') 
in the optical as well as the X-ray region for nickel and iron. For the optical region, we 



show our results in comparison with experiments and the theoretical results of Ref. |T5[ For 
comparison purposes we have not included the Drude term to the diagonal terms in the 
figures for axx{^)- Later, they have been included in the calculation of the magneto-optical 
effects. Fig. 1(a) shows the results of o'l^lu) for nickel. Our theoretical curve seems to be 
in good agreement with the experiments of Ref. BO. The difference between our results and 



those of Ref. ^ is due to a different choice for the lifetime. After repeating the calculation 
with different values of the lifetime parameter, we chose a value of 0.0368 Ry for nickel to give 
best results for MOKE. We also observe the 1 eV shift at 5.4 eV which has been ascribed to 
the failure of the LDA in producing some nickel 3d bands.ill In the case of ujaly{uj) (Fig. 1(b)) 
our results are similar to those of Ref. |15| but the peaks are more pronounced in our case. 
Another feature worth noting is the dip near 5.5 eV that is closer to the observed dip in 
our calculation than theirs. This has a very noticeable effect on the Kerr angle spectra as 
we shall see later. In the case of iron, our theoretical results for crlxi^) (Fig- 1(c)) compare 



far better with experiment than those of nickel and agrees well with the results of Ref. [15. 
Our theoretical results seem to be closer to the experiments of Ref. ^ and agreement in 
general can be taken to be quite good. In the case of u!aly{uj) (Fig. 1(d)) our theoretical 
curve displays a peak at 2.7 eV, that is noticeably higher than the experimentally observed 
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peak. Here, however, the overestimation is surprising since we used an inverse hfe time of 



0.06 Ry for iron, higher than the 0.05 Ry used by Ref. |15. 

Employing the KK transform method outhned earher, we calculated crl^^u) and aly{uj). 
Figure 2(a) shows the result for ual^^u) in the case of nickel while Fig. 2(b) shows the 
one for iron. In both cases, our results are remarkably close to the results of Ref. ITB. In 
the case of nickel they also compare well with experiments.ii In the case of iron the two 
experimental results seem to differ quite widely above 2 eV, thus making comparison with 
theory quite difficult. Figures 2(c) and Fig. 2(d) show our results for ualy^u). For nickel, 
the disagreement at 5 eV is very obvious. In our calculations, this does not seem to have 
affected the polar Kerr angles as seriously as it has the equatorial Kerr effect in the 0-10 eV 
region. Other than this, the results for nickel seem to agree very well with the general 
features of the experiment. For iron, the results of ujal.y{uj) follow the experiment quite well 



and also compare well with the results of Ref. 15. 



B. MOKE and X-MCD 

Using the calculated curves for the elements of the conductivity tensor we evaluated the 
Kerr angles for nickel and iron for the polar Kerr geometry using Eq. (pT]). The resulting 
curves for the optical region are shown in Fig. 3(a) and Fig. 3(b). In the calculation of the 
Kerr angles, we have included the effect of a phenomenological Drude term using values for 
ctd and td from previous experimental results.0 The results on nickel are particularly good 
since the shift of 1 eV is not as noticeable in the Kerr angle spectra as it is in the elements 
of the conductivity tensor. This is in slight contrast to the observations of Ref. |15. 

In a later publicationlll Oppeneer et al have investigated the dependence of MOKE 
spectra on the strength of the spin-orbit coupling. They concluded that the MOKE peaks 
scale linearly with the spin orbit coupling parameter and varying could produce a better 
agreement with the observed MOKE spectra for nickel. No such adjustments for ^ were 
necessary in our results and our theoretical results for nickel seem to agree very well with 
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the experimental results of Ref. ^ throughout the energy range of the data. In the light 
of these observations we may conclude that the spin-orbit coupling strength in our d-bands 
is quite accurate. In the case of iron, our Kerr angles are in excellent agreement with 
experimental results^ and also compare well with the theoretical results of Ref. |T5| and 
Ref. 0. 

As noted before, one of the features of our KK transformation method is that for functions 
that do not decay to zero within the required energy range of interest, accurate momentum 
matrix elements for energies up to two times as much are needed to correctly produce the 
KK transformation. Our momentum matrix elements are sufficiently accurate up to those 
energies since they are calculated using simple analytic expressions resulting from the use 
of a Gaussian basis and hence are free of any numerical approximations. In addition, it is 
evident from our results of Kerr angles that with a straightforward inclusion of spin-orbit 
coupling (for which we again have analytic expressions) in a manner described in the earlier 
section, we are able to very effectively and efficiently account for all the principal features 
of the MOKE spectra. 

It is interesting to see whether the strength of the spin-orbit coupling in our (i-bands 
that has given us good results for MOKE, gives consistent results for the other magneto- 
optical properties as well. This is a very important test since it is now widely acceptedlillii 
that MOKE depends sensitively on the strength of the spin-orbit coupling and the exchange 
splitting. To further investigate whether or not our spin-orbit coupling in the valence bands 
is accurate we decided to calculate the X-ray MOD spectra of iron and nickel at the 2p 
edge using our first principles LOGO method. This serves as a simultaneous check for the 
accuracy of the LDA at X-ray energies. In the X-MCD at the 2p core edge, it is known 
that while the separation between the two peaks originates principally from the spin-orbit 
splitting of the 2p levels, the exact ratio between the two peaks (e.g. approximately -1.6:1 
in the case of nickel) arises out of the spin-orbit splitting of the 3d valence bands. Since X- 
MCD is the difference in absorption of right and left circularly polarized X-rays in the polar 
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geometry, this is nothing but a'^y{uj) evaluated using momentum matrix elements between 
core 2p and valence 3d bands. 

The results of this calculation are shown in Fig. 4(a) for nickel and Fig. 4(b) for iron. 
Along with this, we also computed the al^{uj) which is a measure of the total absorption of 
right and left circularly polarized X-rays. These results are shown in Fig. 5(a) and Fig. 5(b). 
As in any core to band transition the effect of the core hole has to be accounted for, which 
we did by recalculating bands and eigenfunctions with an increased effective Z value and 
then using the new energies and wavefunctions for the final state to calculate the momentum 
matrix elements. The resulting L3 and L2 peaks occur at 856 eV and 874 eV for nickel, which 
agree fairly well with the observed energies of 853 eV and 871 eV respectively. Similarly, 
the calculated L3 and L2 peaks for iron occur at 710 eV and 723 eV, which agree fairly well 
with the observed peaks of 707 eV and 720 eV respectively. However, our main results in 
the X-MCD spectra are the L3-to-L2 ratios which are -1.56:1 in the case of nickel and -1.2:1 
in the case of iron. For nickel this seems to agree very well with the observed ratio of -1.6:1. 
For iron the theoretical ratio is slightly smaller than the observed ratio. However it agrees 
with the ratio for iron determined by other first principle methods based on LDA.ii0 Thus, 
the discrepancy in the L3-to-L2 ratio for iron is not an artifact of our method but may be 
due to the failure of the one-electron band picture. One aspect that must be pointed out 
in the X-MCD spectra of nickel is the missing peaks, often referred to as B and B', 4 eV 
away from the principal peaks on the high energy side. These have been ascribed to many 
body effects in nickel by past author^ and have been reproduced by Jo and Sawatzky in 
a many-body calculation on nickel.0 They have not been reproducible by any one-electron 
band calculation. The results of our band-structure based calculation of MOKE and X-MCD 
are then an indication that, within the LDA framework, our method does give an accurate 
and consistent description of these phenomena. 
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C. The Equatorial Kerr effect and Photoabsorption X-MLD 



The other principal tool for magneto-optical studies on ferromagnetic materials has been 
the equatorial Kerr effect .iHl As mentioned above, the reflection coefficient of an incident 
electromagnetic wave in this geometry depends on the sense of magnetization in the metal. 
When observed at a core-level edge this phenomenon may be termed X-ray magnetic linear 
dichroism. It may be noted that a different kind of absorption X-MLD can be obtained 
by keeping the magnetization constant but rotating the photon polarization vector by 90°. 
This effect has also been calculated previouslyil but is not considered here. 

We have calculated the equatorial Kerr effect for absorption both in the optical region and 
in the soft X-ray region (X-MLD) using Eq. ([121) for several different angles of incidence. 
The results for the optical region are shown Fig. 6(a) and Fig. 6(b) for nickel and iron 
respectively. From Fig. 6(a) it can be seen that as the angle of incidence is changed from 
45° to 80° on nickel the sign of the effect is reversed. After this, from 80° to 88° the peak 
magnitudes progressively increase reaching a peak somewhere close to grazing incidence. 
The effect, of course, disappears at exactly 90°. For nickel, it is important to remember that 
in the evaluation of the conductivity tensor, the theoretical results have always predicted a 
dip at 5 eV where there actually is an experimental rise (see Fig. 1(b) and Fig. 1(d)). Taking 
a clue from this observation, we may predict that in Fig. 6(a), although the structure up 
to 4 eV may compare well with experiment, the peaks (or dips) at about 5 eV may well be 
found to be reversed for every angle. 

The curves of Fig. 6(b) for iron can however be considered to be faithfully reproducing 
experimental results. For iron, we see that as the angle is increased from 80° to 85°, there is 
a reversal of the sign of the effect only for the region upto 5 eV. The effect is most enhanced 
at 88° to 89° and the most observable change is at about 6 eV. Although we did not find any 
experimental results for change in absorption upon reversal of magnetization in the optical 
regime, experimental results for the change in the reflection intensity between magnetized 
and unmagnetized nickel have been reported in the past.ii We calculated this also using 
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Eq. (jl^) by putting Oxy equal to zero to simulate the unmagnetized nickel. Although this 
procedure may be viewed with caution, our results for nickel as shown in Fig. 6(c) agree 
remarkably well with experiment. There again is the characteristic dip at 5 eV which does 
not follow the rise seen in experiment, but this has already been noted to be a consistent 
failure of LDA in nickel. 

For the X-ray region our results for the 2p edge have been shown for nickel and iron 
in Fig. 7(a) and Fig. 7(b) respectively. We clearly see that at the onset of the L2,3 edges 
in both cases the sign of the peaks is the same. There is, however, a sharp overshoot to 
the opposite sign just after the L2 edge. This is more pronounced in the case of nickel, 
although it is unmistakable even in iron. It would be interesting to see if these features are 
actually observed. In the case of iron, the predicted peaks at the L3 and L2 edge seem to be 
of almost equal magnitude. This may be in slight disagreement with experimental results 
since, as mentioned earlier, the prediction of the L3-to-L2 ratios both in the X-MCD and 
total X-ray absorption calculation are slightly underestimated for iron. It may also be noted 
that at the incident angles considered the maximum effect for iron is about 0.8% while that 
for nickel is about 1%. This agrees with the analysis about this type of linear dichroism 
done in the past.@ To the best of our knowledge, such absorption X-MLD experiments at 
the 2p edge for iron and nickel have not been reported although X-MLD in photoemission 
has been subjected to an onslaught of theoretical and experimental analysis, 
results there are characteristically different from our absorption results, since the final states 
in a photoemission experiment are very high above the 3(i bands. As such, the final state 
spin-polarization does not play as significant a role as it does in photoabsorption type mea- 
surements. This is particularly demonstrated in photoemission-XMLDll where the signs of 
the two peaks are opposite each other, in contrast to what is seen in our photoabsorption- 
X-MLD. 
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D. Soft X-ray Faraday rotation 



We finally turn to the Faraday rotation of the plane of polarization of linearly polarized 
X-rays upon transmission through magnetic metals. The main concern here stems from the 
knowledge that the Eq. (|T3|) derived for the Faraday effect is based on the dipole approxi- 
mation. At a first glance, to use it for evaluating the Faraday effect in the soft-X-ray regime 
is to bring the approximation to a questionable limit. This problem has been addressed 
previously^ and it has been demonstrated that the theory of Faraday effects can be ex- 
tended to the X-ray regime by postulating the existence of a free carrier effective dielectric 
medium with the same conductivity tensor as the ferromagnetic metal in question. In fact 
we already did make this extension when we used a'^y to determine X-MCD and cr^^ to de- 
termine the total absorption. Since those results produced reasonably good agreement with 
experiments, we wish to examine its applicabilty further by evaluating the Faraday effect at 
soft X-ray energies. This calculation is a trivial extension to the calculation of X-MCD and 
total absorption since all we need is the KK tranforms of these curves. This, as mentioned 
earlier, can be done very efficiently and accurately with our method. 

The calculated Faraday rotation at the 2p edges of nickel and iron is exhibited in Fig. 8(a) 
and Fig. 8(b). An experiment on soft X-ray Faraday rotation has recently been done at the 
2p edge of iron! with which we compare our results. As can be seen from Fig. 7(b) the 
agreement with experiment is rather good at the L2 edge but is relatively poor at the L3 
edge. This again is related to the underestimated peak in the MCD and total absorption 
spectra at the L3 edge of iron. Nevertheless, given that the dipole approximation is not 
strictly valid at these energies, the agreement with observed angles can be taken to be quite 
good. For the case of nickel, because our X-MCD and total absorption results are in much 
better agreement with experiments, we are more confident of the relative magnitudes of 
angles at the L3 and L2 edges. 
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IV. CONCLUSIONS 



A first principles, self consistent LCGO band calculation has been performed and the 
principal magneto-optical properties have been computed for bulk iron and nickel in this 
paper. In spite of the simplistic way in which spin-orbit coupling has been included, our 
results have agreed very well with previous first principles calculations based on LDA of 
MOKE and MCD and experimental observations. In particular our MOKE results for nickel 
have produced very good agreement with experiments, despite the lack of good description 
of some of the LDA bands in nickel. Discrepancies originating from the failure of the LDA 
have however crept up in the equatorial Kerr effect in the optical region around 5 eV. Our 
results for photoabsorption X-MLD could not be compared with experiments but judging 
from the agreement of the X-MCD and soft X-ray rotation results with experiments, we 
conclude that our results for nickel may be considered to be accurate, but for iron the actual 
L3 peak magnitude may be much higher as compared to the L2 peak. This conclusion is also 
supported by our Faraday effect results on iron. It is also noteworthy that the equatorial 
Kerr effect is larger in the optical region than in the soft-X-ray region. 

We have also demonstrated that fast, efficient and accurate KK transforms can in fact 
give sufficiently satisfactory results for sensitive magneto-optical effects of bulk metals, which 
have been previously obtained by computationally intensive methods. This speed and ef- 
ficiency is particularly important if one is to use such calculations for quickly obtainable 
results for technological applications. This analysis has also confirmed that expressions for 
the Faraday effect for the optical region may be extended to the soft X-ray region yielding 
satisfactory results. The principal feature of this work is however that a single calculation 
of all the components of the conductivity tensor in both the optical as well as the X-ray 
region has yielded consistent results for a host of different magneto-optical effects. 
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FIGURES 

FIG. 1. Elements of the conductivity tensor for nickel and iron. In (a) and (b), solid line with 



an inverse lifetime of 0.0368 Ry is our result. Dashed line is the result of the calculation of Ref. 15 
with an inverse lifetime of 0.04 Ry. In (c) and (d) solid line is our result for an inverse lifetime of 
0.06 Ry. Dashed line is the result of the calculation of Ref. ^ for an inverse lifetime of 0.05 Ry. 



(a) a^.^ for nickel. Circles are the experimental results of Ref. 50. Squares are the experimental 



results of Ref. ^ (b) Loa'^y for nickel. Circles are the experimental results of Ref. (c) cr^^ 



for iron. Circles are the experimental results of Ref. 52. Squares are the experimental results of 



Ref. 55. (d) uja'^y for iron. Circles are the experimental results of Ref. 64 



2 

XX 



FIG. 2. Dispersive components of the conductivity tensor for nickel and iron, (a) ua. 
for nickel. Circles are the experimental results of Ref. (b) wo"^.^. for iron. Circles are the 
experimental results of Ref. Squares are the experimental results of Ref. (c) ujal.y for nickel. 
Circles are the experimental results of Ref. (d) u>al.y for iron. Circles are the experimental 



results of Ref. 64 



FIG. 3. Polar Kerr rotation. In both graphs, solid line is our result. Dashed line is the result 
of the calculation of Ref. |l^. (a) For nickel. Circles are the experimental results of Ref. Squares 
are the experimental results of Ref. ^ (b) For iron. Circles are the experimental results of Ref. |5^. 

FIG. 4. X-MCD at the 2p edges for (a) nickel and (b) iron. A Gaussian broadening of 0.2 eV 
has been added in both cases. 

FIG. 5. Total absorption at the 2p edges for (a) nickel and (b) iron. A Gaussian broadening 
of 0.2 eV has been added in both cases. 
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FIG. 6. Equatorial Kerr effect. For (a) and (b) AT is the cliange in absorption intensity upon 
reversal in magnetization. T is the average absorption of the two directions. Plots are correct upto 
a negative sign, (a) In absorption for nickel and (b) In absorption for iron. = 45° Solid line. 
e = 80° Dotted hne. 9 = 85° Dashed line. 9 = 88° Dash-dotted hne. (c) 61 = 45°. A/ is the change 
in reflection intensity upon magnetization. / is the reflection intensity from unmagnetized nickel. 
Solid line is our theoretical result. Circles are the experimental results of Ref. 

FIG. 7. Photoabsorption X-MLD at the 2p edges. AT is the change in absorption intensity 
upon reversal in magnetization. T is the average absorption of the two directions. Plots are correct 
upto a negative sign, (a) for nickel. 9 = 85° Solid line. 9 = 87° Dotted line. 9 = 88° Dashed line, 
(b) for iron. 9 = 80° Solid line. 9 = 85° Dotted hne. 9 = 87° Dashed hne. 

FIG. 8. Soft X-ray Faraday rotation at the 2p edges. In both figures solid line is our theoretical 
results, (a) for nickel, (b) for iron, (i = 80 nm. Circles are the experimental results of Ref. |^. 
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